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We determine the zero-bias anomaly of the conductance 
of tunnel junctions by an approach unifying the conventional 
Coulomb blockade theory for ultrasmall junctions with the 
diffusive anomalies in disordered conductors. Both, electron- 
electron interactions within the electrodes and electron-hole 
interactions between the electrodes are taken into account 
non-perturbatively. Explicit results are given for one- and 
two-dimensional junctions and the crossover to ultrasmall 
junctions is discussed. 
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Coulomb blockade of tunneling has been investi- 
gated extensively in the last decade. Most of these stud- 
ies were focussed on ultrasmall tunnel junctions with 
small capacitance C. Then, the relevant energy scale 
is the charging energy Ec = e^/2C. In these effectively 
zero-dimensional junctions the charge of a tunneling elec- 
tron is distributed over the junction area fast compared 
to the time scale h/eV associated with the applied volt- 
age V, and Coulomb blockade effects are observable for 
single junctions only because of the slow charge relax- 
ation via the environmental impedance Z(io) |2|. 

The situation changes for tunnel junctions with large 
lateral extensions where the time scale of propagation of 
the electromagnetic field over the system size becomes 
relevant. Then screening and relaxation processes within 
the electrodes need to be considered explicitly. Particu- 
larly, in junctions with disordered low-dimensional elec- 
trodes, the field propagates diffusively, and a suppression 
of the tunneling conductance arises from the reduction of 
the electron density of states by electron-electron inter- 
actions [H]. 

The pioneering papers on the diffusive anomalies in dis- 
ordered conductors are based on perturbation theory 
in the Coulomb interaction and fail at very low energy 
scales. As is well-known from the conventional theory of 
Coulomb blockade in effectively zero-dimensional junc- 
tions a non-perturbative approach is needed for a 
complete treatment of the Zero-Bias Anomaly (ZBA) of 
the conductance. Furthermore, apart from the interac- 
tion within the electrodes, the Coulomb interaction be- 
tween the electrodes needs to be taken into account. Only 
in limiting cases, e.g., for a strongly disordered electrode 
screened by a bulky electrode, can the mutual interaction 
be incorporated into a screened Coulomb potential. 

Here, we present a microscopic approach treating 
all low energy excitations of the Coulomb field non- 
perturbatively and obtain explicit results for the current- 



voltage relation of zero-, one-, and two-dimensional 
junctions at finite temperature. Specifically, we consider 
a tunnel junction consisting of two disordered metallic 
electrodes separated by a thin insulating barrier with a 
chemical potential difference /ii — = e^- The tun- 
neling current / can be calculated from the tunneling 
Hamiltonian 

i/^ = ^Tfc,clc, + h.c., (1) 

where k and q label the momenta of electrons in the right 
and left electrode, respectively ||]. To leading order in 
Ht the current reads 

I{V) = -2elm[X'''\eV)]. (2) 

Here is the Fourier transform of 

fc.g k' ,q' 

where we put h — I, and the time evolution of the Fermi 
operators arises from the Hamiltonian in the absence of 
tunneling, but in the presence of disorder and interac- 
tions. 

To evaluate the expectation value with four Fermi op- 
erators in Eq. (||) , we employ the Keldysh cr-model tech- 
nique 11,0] leading to the result in Eq. ([l^ ) below. We 
start from the partition function 



Z =JV VipiVipiVip2'Dil;2'D(l) e 



iS[tl>l,tl>l,tl>2,i>2,'P] 



(4) 



where A'' is a normalization constant. The action is given 
by 



J2 f dt [ dr4,,{r,t) 
. Jc JVi 



i— + —A 
dt 2m 



-Wj{r) + e(l3{r, t) 



iPjir, t) 



-— [ dt I dr 



V0(r, t) 



(5) 



where C is the Keldysh contour, Wj{r) are Gaussian dis- 
tributed disorder potentials, 4> is the electric potential, 
and e is the electric charge. While the fermionic terms 
are integrated over the volume of the electrodes only, the 
last term is integrated over entire space, which leads to 
electrostatic coupling between the electrodes. For sim- 
plicity we disregard polarization effects in the tunnel bar- 
rier, that may easily be incorporated in the capacitance 
Co introduced below. 
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Assuming that the impurity potentials in the two elec- 
trodes are statistically independent, the Gaussian inte- 
grations over the realizations of Wj{r) lead to quartic 
interactions of the fermion fields within each electrode. 
Within the cr-model approach these interactions can be 
decoupled by introducing Hubbard-Stratonovich matrix 
fields Qj 1^ . Subsequently, the fermionic fields can be in- 
tegrated out in the usual way yielding a representation of 
the impurity averaged partition function as a functional 
integral over the electric potential 4> and the matrix fields 
Qj. Since for given field 4> all of these transformations 
can be done independently for each electrode, the result- 
ing action S'[0, Qi, Q2] can immediately be inferred from 
earlier work 

To proceed we follow Kamenev and Andreev Q and 
evaluate the integrals over the matrix fields Qj using the 
saddle point approximation. The saddle point solution 
can be obtained analytically for spatially uniform fields 
(j){t) only. Provided (f> is essentially constant on scales of 
the size of the mean free path I, an effective action was 
derived in which incorporates dynamic screening in 
the random phase approximation (RPA) and the diffuson 
vertex correction within each electrode. 

One is left with a Gaussian integral over the Coulomb 
field with a purely electromagnetic action which deter- 
mines the field propagators. Since these depend on the 
geometry of the electrodes, we first restrict ourselves to 
the case of thin films. Then, for low energy excitations 
with wave vectors q ^ a,~^, where a is the thickness of 
the films, we have effectively a two-dimensional prob- 
lem, where the Coulomb field depends on two spacial 
coordinates in each electrode (see also the discussion be- 
low). Although only the large distance behavior of the 
Coulomb interaction matters, we have to keep the bar- 
rier thickness A finite since dipolar interactions arise . 
With the sources restricted to two-dimensional films at 
Zj = ±A/2, the bare Coulomb interactions reads 



x' ,y-y') 



dzdz 



, e'^S{z - Zi)S{z' - Zj) 



The resulting equation for the matrix f7''°'(q, w) of the 
Fourier transformed field propagator then takes the form 
of the diffusive RPA 



where 



(6) 



(7) 



Here u{q) and v(q) are the two-dimensional bare 
Coulomb interactions inside and between the electrodes, 
respectively. Neglecting polarization effects in the bar- 
rier, we obtain u{q) = 27re^/g and v{q) 
The polarization function 



2Tre^e"J^/q. 



Piq,io) = 



_Dj£_ 



V2 







(8) 



contains the bare electron densities Vj and the electron 
diffusion constants Dj of the two electrodes. 

Supplementing the partition function Z with source 
fields that couple to the fermion operators of each elec- 
trode, it may be employed as a generating functional to 
evaluate impurity averaged correlation functions like 



k,q k\q' 



kq k' q' 



4Wc,(t)4c,,). (9) 



The pre-exponential factors resulting from functional 
derivatives with respect to the sources are replaced by 
their saddle point approximation. In terms of an average 
tunnel matrix element T, we find 



X>{t) = -iyiU2\T\^ j de J de' e^(^-^')* 
xn{e)[l ~ n{e')]e''^'\ 
where the interaction effects are incorporated in 



Jit) 



duj e 



27r 1 



1 



Iniy(w) 



Here (3 is the inverse temperature and 

Uij {q, 



{D^q^ — iuj){Djq'^ — iuo) 



(10) 



(11) 



(12) 



The factors [Djq^ ~ iuj)^^ come from the vertex correc- 
tion terms in the action and are seen to obey the classical 
diffusion equation. When the two-dimensional result for 
y(a>) is inserted into Eq. (pT|), the cj-integral needs to be 
cut off at frequencies of order \uj\ « I/tq ~ a^/D where 
the crossover from two- to three-dimensional behavior 
occurs. 

Now, combining Eqs. (||) and (|l0|), the expression for 
the current is found to read 



liV) 



Gq 



de eP{eV -e)- 



-/3eV 



1 - e-/3' 



(13) 



where Go = 47re^:/iz/2|rp is the bare tunneling conduc- 
tance, and where we have introduced the spectral density 



1 r 

P{e) = -Re / 
7^ Jo 



(14) 



The expression ( [l^ ) constitutes the central result of this 
work. It is formally identical to the conventional expres- 
sion for the current -voltage relation of ultrasmall tunnel 
junctions (zero-dimensional case) and determines the 
non-perturbative effect of Coulomb interactions on the 
I-V-cuTve of spatially extended disordered tunnel junc- 
tions. The very same form can also be derived in the case 
of effectively one-dimensional wires (see below) . 
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To make contact with previous work, we first con- 
sider the (somewhat unreahstic) case where Coulomb in- 
teractions between the electrodes are disregarded, i.e., 
Ui2 = U21 = 0. Then J{t) = Ji(t) + J2{t) where J -jit) 
contains the contribution from Ujj only. Replacing J(t) 
by Jj{t) in Eq. (p^, we may define spectral densities 
Pj{e) for each electrode. The current (13) can then be 
written in the familiar form 



/ 47re|r| 



de vi{€)v2{e - eV) [n{e - eV) - n{e)] (15) 



with the densities of states 

/•OC 



de' 



1 



1 



(16) 



This relation gives a non-perturbative result for the effec- 
tive tunneling density of states of an electrode in presence 
of Coulomb interactions expressed in terms of the spec- 
tral density Pj (e) familiar from Coulomb blockade theory. 
When -Pj (e) is replaced by its perturbative approxima- 
tion, the seminal result by Altshuler, Aronov, and Lee 
is recovered. 

In particular, for a two-dimensional film at zero tem- 
perature the density of states ( [l^ ) reads 



z/(e) = z/exp 



1 

47r5 



ln|eTo|ln 



(17) 



where we have suppressed the index j. Further, g = 
2'nvD is the conductance in units of e^/27r and k = 
is the inverse screening length in two dimensions. The 
non-perturbative result (|l^) has been obtained previ- 
ously by Kamenev and Andreev [Q. 

We now return to the full problem and determine the 
non-perturbative I-V-cxvcve for two-dimensional inter- 
acting electrodes at zero temperature. Then, the four 
terms in Eq. (12) cannot be split into contributions of 
each electrode, and the spectral function P(e) character- 
izes the coupled system. In the parameter range relevant 
for tunneling experiments u{q) — v(q) « 27re^A, and we 
obtain from Eq. (f2|) 



F(^) 



e 



EE 



1 1 



^ {J)*q^ — iLu){Djq'^ — iuj) 



(18) 



Here we have introduced the field diffusion constant D* ~ 
{2SiS2+SiD2 + S2Di)/{Si+52), with Sj = DjKjA and the 
numerical factors Xj = {2Si + Di — Dj)/2{Si+62) {i j)- 
While D and k are properties of the electrodes, A may 
also be expressed in terms of the capacitance per unit 
area Co = l/47rA of the junction. With the result ( p^ ) 
we readily obtain from Eq. ([l^) for the zero temperature 
differential conductance G{V) ~ dl/dV at voltages V <^ 
Vo 



-27/9 



V 



GjV) ^ 
Go T{l + 2/g)\Vo 



2/s 



(19) 



where 7 = 0.577... is Euler's constant, g ~ 
^CoD* a, ln(l/C,)/(l - with = D,/D* 

is a dimensionless parameter, and Vo — 27r/eTo « 
27rD/ea2. 

When g ^ 1 we recover the perturbative result 



G{V)^Go[l + 2/g\n{V/Vo)], 



(20) 



Hence in the non-perturbative approach the logarithmic 
corrections of perturbation theory are exponentiated to a 
power-law dependence on V. In the limiting case where 
one of the electrodes is bulky, i.e. ki S> K2, this result can 
be rewritten by means of Eq. ( p^ ) in terms of an effective 
tunneling density of states of the other electrode 



-27/9 



m + 2/g) 



e 



2/9 



(21) 



This gives a non-perturbative generalization of the result 
obtained by Altshuler, Aronov, and Zyuzin Q. 

We now turn to one-dimensional contacts. Then the 
result (|l8|) for Y{uj) remains valid provided the bare in- 
teractions u{q) and v{q) in Eq. (Q) are replaced by u{q) — 
2e^ln(l/ga) and v{q) = 2e'^ \n{l / qA) . Further the mo- 
mentum sum is one-dimensional, and Co — [41n(A/a)]^^ 
becomes the capacitance per unit length. The zero tem- 
perature differential conductance now takes the form 



GjV) 
Go 



= 1 - erf A 



with Vq 



8^E,A,/(i+er)F. 



(22) 



This result is 



formally equivalent to the conductance of an ultrasmall 
junction biased via a RC transmission line | [ic| ] and is in 
quantitative agreement with a recent experimental study 
of long tunnel junctions [ pT| . 

If one of the electrodes is bulky, we again may use 
Eq. (^6|) to obtain an effective tunneling density of states 
of a one-dimensional electrode 



u{e) = 1/0(1 — erf- 



eVo 



(23) 



While for large energies we recover the perturbative re- 
sult vje) — fo(l — 2-\/eVo/7re), the non-perturbative re- 
sult (E3) does not diverge at low energies but approaches 
an exponential suppression of the densities of states near 
the Fermi surface ^{e) ~ vq\J tj-neVo exp(— eVb/e). 

For the realistic case of a system with finite size the 
diffusive spreading of the transfered charge reaches the 
boundaries of the electrode for long times, and then the 
charge relaxes via the external circuit characterized by 
an impedance Zext(w). This modifies the ZBA at very 
low voltages. For simplicity we assume that the external 
charge relaxation, which is the only process relevant for 
ultrasmall junctions O, is slow and couples only to the 



3 



q ~ component of the field 
modified to read 



12|. Then Eq. (pt) is 



A, 



9#0 j- 



iijj){Djq'^ 



iuj) 



(24) 



For one-dimensional electrodes of length L the discrete q 
values are of the form (27r/i)n, n integer, and C = CqL 
is the total capacitance of the contact, while for two- 
dimensional quadratic electrodes C = CqL"^ . 

For three-dimensional electrodes of finite size there 
are several energy scales in Ylu): The charging energy 
Ec = e^/2C and for each spatial dimension i of the 
electrodes a Thouless energy ii^Th = D/i"^ as well as a 
field Thouless energy E* — D* /P. Here I stands for the 
length L, width w, and thickness a, respectively. The ef- 
fective dimension of the electrodes is determined by E* . 
Whenever eV exceeds one of the scales i?! [H — L, w, a) 
the effective dimension increases by one. 

As an example we consider a symmetric one- 
dimensional contact of length L with an ohmic external 
impedance Z^^t = R- Then, from Eq. (124) we obtain 



(25) 



l + a2(tji?KC)2 



where luvYoo{u)) = /[2\/2D*Cq{\ + ^1/2)^3/2] de- 
scribes an infinitely long contact, and the function 
f{z) = {g{z) - e'^9{z/0)/{l - e'^) with g[z) ^ 
Ini [— \/2i cot (Viz/2)] incorporates all size effects. The 
last term in Eq. (E5|) describes the environmental 
impedance characterized by the dimensionless parame- 
ter a = R/Rk with Rk = h/e^. 

In Fig. |l| we present results for typical values of metal- 
lic quasi one-dimensional tunnel junctions (^ = 
4 X 10-'', E*/eVo = 100, Ec/eVf^ = 25). The full 
line corresponds to a finite length junction shunted by 
a rather large environmental impedance (a = 1) to dis- 
play the crossover between zero-dimensional and one- 
dimensional behavior clearly. The dotted line gives the 
conventional environmental Coulomb blockade of a zero- 
dimensional junction which determines the behavior for 
voltages below Ec/e (see inset). On the other hand, for 
large voltages above E* /e the result ( p2| ) for a long one- 
dimensional junction (dashed line) is approached. The 
chain dotted line results from a perturbative calculation 
of the diffusive anomaly . 

In summary, we have studied tunneling in large junc- 
tions with diffusive electron motion. Treating the 
Coulomb interaction non-perturbatively, the /-1^-curve 
was written in a form familiar from Coulomb blockade 
theory for ultrasmall junctions. Interaction effects are 
again incorporated in a spectral function P(e). However, 
while for ultrasmall junctions P(e) depends only on the 



effective impedance [Z~^^{uj) — iujC]~'^ , for large junctions 
this quantity is replaced by —iu!Y{u!)/e'^ which incorpo- 
rates also the diffusive anomaly of disordered metals. 
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FIG. 1. Differential conductance G(y) of a quasi 
one-dimensional tunnel junction of finite length with an ex- 
ternal impedance at T = (solid line), see text for details. 
This is compared with the results of a zero-dimensional (dot- 
ted line) and a long one-dimensional junction (dashed line). 
The chain dotted line corresponds to the perturbative result. 



[1] Single Charge Tunneling, Vol. 294 of NATO ASI Series 
B, edited by H. Grabert and M. H. Devoret (Plenum, 
New York, 1992). 

[2] M. H. Devoret, D. Esteve, H. Grabert, G.-L. Ingold, H. 
Pothier, and C. Urbina, Phys. Rev. Lett. 64, 1824 (1990); 
S. M. Girvin, L. I. Glazman, M. Jonson, D. R. Penn, and 
M. D. Stiles, Phys. Rev. Lett. 64, 3183 (1990). 

[3] B. L. Altshuler, A. G. Aronov, Sov. Phys. JETP 50, 968 
(1979). 

[4] B. L. Altshuler, A. G. Aronov, P. A. Lee, Phys. Rev. 

Lett. 44, 1288 (1980). 
[5] We restrict ourselves to spinless fermions here but include 

spin degeneracy by appropriate factors of 2. 
[6] D. Belitz, T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 

(1994). 

[7] A. Kamenev, A. Andreev, Phys. Rev. B 60, 2218 (1999). 
[8] A is actually the dipole moment per unit charge which 
due to the finite film thickness is somewhat larger than 
the actual barrier thickness. 
[9] B. L. Altshuler, A. G. Aronov, A. Yu. Zyuzin, Sov. Phys. 
JETP 59, 415 (1984). 
[10] G.-L. Ingold and Y. V. Nazarov, in Ref. 1, p. 21. 
[11] F. Pierre, H. Pothier, P. Joyez, N. O. Birge, D. Esteve, 

and M. H. Devoret, Phys. Rev. Lett. 86, 1590 (2001). 
[12] A. Kamenev, Y. Gefen, Phys. Rev. B 54, 5428 (1996). 



4 



